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Abstract — The aim of this paper is to determine quantum 
master and filter equations for systems coupled to continuous- 
mode single photon fields. The system and field are de- 
scribed using a quantum stochastic unitary model, where the 
continuous-mode single photon state for the field is determined 
by a wavepacket pulse shape. The master equation is derived 
from this model and is given in terms of a system of coupled 
equations. The output field carries information about the system 
from the scattered photon, and is continuously monitored. The 
quantum filter is determined with the aid of an embedding 
of the system into a larger system, and is given by a system 
of coupled stochastic differential equations. An example is 
provided to illustrate the main results. 

I. Introduction 

In recent years single photon states of light have be- 
come increasingly important due to applications in quantum 
technology, in particular, quantum computing and quantum 
information systems, [20], [23], [18], [14], [26]. For instance, 
the light may interact with a system, say an atom, quantum 
dot, or cavity, and this system may be used as a quantum 
memory, [20], or to control the pulse shape of the single 
photon state [23]. Note that in practice one can consider 
different types of 'single photon states', including the single 
photon state of a single mode of light confined inside an 
optical cavity [20], or a single photon state superposed over 
a continuum of modes of a travelling field (i.e., colloquially, 
a "flying" single photon state) referred to as a continuous- 
mode or multimode single photon state [19], [22], [23]. In 
this paper we will be interested exclusively with the latter 
kind of single photon state and therefore from this point 
on when we say 'single photon' state we specifically mean 
the continuous-mode single photon state of a travelling field. 
When light interacts with a quantum system, information 
about the system is contained in the scattered light. This 
information may be useful for monitoring the behavior of the 
system, or for controlling it. The topic of this paper concerns 
the extraction of information from the scattered light when 
the incoming light is placed in a single photon state, denoted 
1^), as illustrated in Figure [T] 

The problem of extracting information from continuous 
measurement of the scattered light is a problem of quantum 
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Fig. 1. System G coupled to a field in a continuous-mode single photon 
state 15) defined by the wavepacket ^. The output field is continuously 
monitored by homodyne detection (assumed perfect) to produce a classical 
measurement signal Y{t). 



filtering, [4], [5], [6], [10], [27], [3], [8], [28]. The current 
state of the art for quantum filtering considers incoming 
light in a vacuum or other Gaussian state, with quadrature 
or counting measurements. Single photon states of light are 
highly non-classical, and are fundamentally different from 
Gaussian states. At present, to the best of our knowledge, 
there are no filtering results for systems driven by single 
photon fields. In view of the increasing importance of single 
photon states of light, the purpose of this paper is to solve 
a quantum filtering problem for systems driven by single 
photon fields. As a by-product, we obtain a master equation 
for the system. A significant feature of the master equation 
and the quantum filter is that they are both given by a system 
of coupled equations, not a single equation as in the vacuum 
case. This reflects the non-Markovian character of systems 
driven by single photon fields. 

The paper is organized as follows. The filtering problem 
to be solved is formulated in Section |ll] The master equation 
is derived in Section|lII]using the model presented in Section 
Hn This leads naturally to Section |IV] where the system 
is embedded in a larger Markovian model, using a signal 
generator model. The single photon filter is presented in 
Section |V] and an example is discussed in Section |VI] 

Notation: We use the standard Dirac notation \ip) to 
denote state vectors (vectors in a Hilbert space) [21], [1]. 
The superscript * indicates Hilbert space adjoint or complex 
conjugate. The inner product of state vectors |-0i) and \ip2) is 
denoted {ijji\ip2)- The expected value of an operator X when 
the system is in state 1-0) is denoted E^[X] = {ip\X\ip). For 
operators A and B we write (A, B) = tr[A*B]. 

II. Problem Formulation 

We consider a quantum system G coupled to a quantum 
field Bin, as shown in Figure [T] The interaction of Bin with 
G produces the output field Bout- The input field is placed 
in a single photon state, denoted using Dirac's notation 



as 1^), where ^ is a complex valued function such that 
/o IC(-5)P'^s = 1- As illustrated in Figure [T] the wavepacket 
1^) interacts with the quantum system G, and the results of 
this interaction provide information about the system that 
may be obtained through continuous measurement of an 
observable Y{t) of the output field Bout{t)- The filtering 
problem of interest in this paper is to determine the condi- 
tional state from which estimates X{t) of system operators 
X may be determined at time t based on knowledge of the 
observables {Y{s), < s < t}. 

In what follows the system G is assumed to be defined on 
a Hilbert space ^5, with an initial state denoted \rj) e S)s- 
The input field Bin is described in terms of annihilation 
B{^) and creation B*{£) operators defined on a symmetric 
(Boson) Fock space ^, [24, Chapter II], [8, Section 4]. 
The continuous-mode single photon state is defined on the 
symmetric Fock space by [19, sec. 6.3], [22, Section 14.2], 
[23, eq. (9)] 

IO = s*(OI^>, (1) 

where |0) is the vacuum state of the field. Expression ([T]i 
says that the single photon wavepacket is created from the 
vacuum using the field operator B*{£_). 

The Hilbert space for the composite system is 

where here we have exhibited the continuous temporal tensor 
product decomposition of the Fock space d = St]'^ S{t into 
past and future components, which is of basic importance 
in what follows. We use the notation E to denote quantum 
expectation, usually with a subscript to denote the state being 
used. In particular, we write 

En[X ® F] = (77^|(X ® F)\r,0 = {v\X\v){(\F\0 (2) 

for the expectation with respect to the product state \rj^), 
where the field is in the single photon state. Here and in what 
follows X is a bounded system operator acting on Sjs, and F 
is a field operator acting on the Fock space ^. Similarly, we 
may define the expectation when the field is in the vacuum 
state, 

Eoo[X®F] = M(X®F)|7?,/)) = (7,|X|77)(0|F|(A). (3) 

We will also have need for the cross-expectations 

Ew[X(g,F]^{r^^\{X(g,F)\ri(j)), 



Eoi[X®F] = (#|(X®F)|r;0. 



(4) 



A crucial difference between the single photon state and 
the vacuum state is that the later state factorizes 10) — \(j)t])<E} 
|(/)(() with respect to the temporal factorization ^ = g'j] ®5^(t 
of the Fock space, with \<f>t]) G St] and \(j)(^t) G T?(t> while 
the former does not. Rather, we have 



10 = B*{om = 1^,]) ® !(/.(,) + 10,]) ® |C(,), (5) 



where 



and i?r(0 = Bi^xm), B+iO = S(Cx(t,oo)), 5(0 = 
Bf. (0 + B^{£_). Here, X[o,t] is the indicator function for the 
time interval [0,i]. Note that while |0 has unit norm, we 
have 



I6]>|P= f\m\'ds, and 11 10* 
Jo 



\m\'ds. 

(7) 

A consequence of the additive decomposition (|5]l is the 
following. Let K{t) be a bounded operator acting on the 
full Hilbert space io that is adapted, i.e. K{t) acts trivially 
on 5^((, the field in the future. Then the expectation with 
respect to the single photon field may be expressed in terms 
of the vacuum state as follows; 

En[Kit)] = Eoo[BriOl<{t)Br{0+wit)K{t)](8) 

where w{t) = J^ \^{s)\'^ds. 

The dynamics of the system will be described using the 
quantum stochastic calculus, [17], [11], [24], [12], [8]. Quan- 
tum stochastic integrals are defined in terms of fundamental 
field operators B{t), B*{t) and A(i), [24, Chapter II], [8, 
Section 4]1^ The non-zero Ito products for the field operators 
are 

dB{t)dB*{t) = dt, dB{t)dA{t) ^ dB{t), 
dk(t)dK{t) = dA{t), dA{t)dB*{t) = dB*{t). (9) 

The dynamics of the composite system is described by a 
unitary U{t) solving the Schrodinger equation, or quantum 
stochastic differential equation (QSDE), 

dU{t) ={{8- I)dA{t) + LdB*{t) - L*SdB{t) 

^{]-L*L + iH)dt}U{t), (10) 

with initial condition C/(0) = /. Here, iJ is a fixed self- 
adjoint operator representing the free Hamiltonian of the 
system, and L and S are system operators determining 
the coupling of the system to the field, with S unitary. In 
this paper, for simplicity we assume that the parameters 
G = {S, L, H) are bounded operators on the system Hilbert 
space S^s- 

A system operator X at time t is given in the Heisenberg 
picture by X{t) ^ jt{X) = U (t)* {X (g) I)U (t) and it follows 
from the quantum Ito calculus that 



djtiX) 



where 



and 



jt{S*XS - X)dA{t) + jt{S*[X,L])dB{ty 
+3t{[L* ,X]S)dB{t) + jt{g{X))dt, (11) 



g{x) = CL(x)-i[x,H], 



(12) 



10]) = Sr*(e)l0t]), and 10*) = i?*+*(OI</'(t), (6) 



CL{X) = h*[X,L] + ^[L*,X]L. (13) 

' In terms of annihilation and creation wliite noise operators b{t),b*(t) 
that satisfy singular commutation relations [b{s) , b* (t)] = 5(t — s), the 
fundamental field operators are given by B{t) = J„ b{s)ds, B*{t) = 
Jq b*{s)ds, and A(i) = /„ b*(s)b{s)ds. Also, we may write -B(5) = 



and The map X i-^ GiX) is known as the Lindblad 
generator, while the quartet of maps X n- Q{X), S*XS — 
X, S*[X,L], [L*,X]S are known as Evans-Hudson maps. 

The output field is defined by Bout{t) = U*{t)B{t)U{t)E 
In this paper we consider the output field observable Y{t) 
defined by 

Yit) = U*{t)Zit)Uit), (14) 



where 



Z{t) = B{t) + B*it), 



(15) 



is a quadrature observable of the input field. Note that 
both Z{t) and Y{t) are self-adjoint and self-commutative: 
[Z{t),Z{s)] = and [Y{t),Y{s)] = 0. We write ^t and 
^ for the subspaces of commuting operators generated by 
the observables Z{s), Y{s), < s < t, respectively^ 
They are related by the unitary rotation ^ = U*{t)3ftU{t). 
Physically, Y{t) may represent the integrated photocurrent 
arising in an idealized (perfect) homodyne photodetection 
scheme, as in Figure [T] For further information on homodyne 
detection, we refer the reader to the literature; for example, 
[2], [3], [28], [16]. In particular, [16] considers pulsed ho- 
modyne detection for fields in a continuous-mode 7i-photon 
state, which includes the single photon state as a special case. 
The primary goal of this paper is to determine the quantum 
filter for the quantum conditional expectation (see, e.g. [8, 
Definition 3.13]) 



X(f)=Eu[X(i)|^]. 



(16) 



This conditional expectation is well defined, since X{t) 
commutes with the subspace ^ (non-demolition condition). 
The conditional estimate X{t) is affiliated to ^ (written in 
abbreviated fashion as X{t) E ^) and is characterized by 
the requirement that 



for all K e^f 



in[Xit)K]^En[X{t)K] 



III. Master Equation 



(17) 



Before deriving the quantum filter, we work out dynamical 
equations for the unconditioned single photon expectation. 
Such equations are often called master equations and are 
of fundamental importance, and arise in Markovian models 
of open quantum systems, [12], [24], [10], [28]. Master 
equations are analogous to the Fokker-Plank equations for 
classical diffusion processes. Note that the master equations 
for systems driven by a single photon field have previously 
been derived by other means in [13], although we only 
became aware of this after this work was completed. 

When the field is in the vacuum state |(/)), the joint 
system-field state evolves according to U{t)\r](l)), and the 
system density operator p"*^ (f ) is defined by (p"° (t) , X) = 



-Recall B{t) = Bi„(t) is the input field. 

^^t and {% are commutative von Neumann algebras. They are also 
filtrations, e.g. J^n C Si^ whenever ti < t2. 



{'i-l(j),U*{t){X (g) I)U{t)-q(j)) = Eoo[X(<)]. It is well-known 
[24], [17], [12] that p'^^{t) satisfies the master equation 



pOO(t)^g*(p"0(i)), 



(18) 



where 



g*{p) = LpL* - ^pL*L - ^L*Lp + i[p, H]. (19) 

The master equation ( fTSI l is readily determined from the 
Heisenberg evolution (fTTT i by taking expectations with re- 
spect to the vacuum state and appropriately collecting terms. 
Note that the unitary operator S appearing in the Schrodinger 
equation ( fTOl i does not appear in the master equation ( fTsT l. 

Now suppose that the field is in a single photon state 
1^), in which case the density operator p^^{t) is defined 
by (p"(t),X) = {rj^,U*{t){X®I)Uit)f^O =Eii[X(t)], 
which involves expectation with respect to the single photon 
field. Using equation (fTTT i and the relations 

dB{t)\o = atM), dA{tm = mdB*{tM), no) 

we calculate that 

lE,^[X{t)] = jEn[Xit)] (21) 

= En[g{Xm 
+Eo,[S*it)[Xit),Litm*{t) 
+EM[L*{t),X{t)]Sit)m 
+Eoo[{S*{t)X{t)S{t)^X{t)mit)\^. 

Notice that the right hand side of (|2T]) includes a vacuum 
expectation, as well as cross terms involving single photon 
and vacuum states. The system driven by the single photon 
field is not Markovian, in contrast to the vacuum case. 
In view of this, we define 



fi','iX) = En[X(t)], fil^iX) = Eio[X(t)] 
p'i\X)=Eoi[X{t)] M™(X)=Eoo[X(i)]. 



(22) 



Consequently, the master equation in Heisenberg form for 
the system when the field is in the single photon state |^) is 
given by the system of equations 

pI\X) = f,l\g{X))+f,<i\S*[X,L])C{t) (23) 

+pl'{[L*,x]s)m 

+p'i'{s*xs~xMtr, 
pl'ix) = ^,l'>igix)) + ^,°,'>is*[x,L])C{t), (24) 

A?^(x) = f,',\gix)) + p°,%[L*,x]s)m, (25) 

A°°(X) = priGiX)). (26) 

The initial conditions are 

/.Ji(X) = p",'{X) = (r,,Xr;), ^,1"{X) = p°HX) - 0. 

(27) 
In order to obtain a Schrodinger form of the master 
equations, we define (generalized) density operators p'^it) 
by 

{p>\t),X)^p.l\x). (28) 



The operators p-' [t) enjoy the symmetry properties 

p°°*(t) = p°°(i), p"^*(t) = pi°(t), p"*(i) = p'\t). 

(29) 
The master equation in Schrodinger form for the system 
when the field is in the single photon state |^) is given by 
the system of equations 

+ [L,p'\t)S*]C{t) 

+{Sp°\t)S* - p^\t)mt)\\ (30) 

p'\t) = g*(pi"(f)) + [^pOO(f),i*]C(t), (31) 

fP\t) = g*{p'>\t)) + [L,p^\t)S*]C{t), (32) 



p°°(i) = a*(p""W). 
The initial conditions are 



(33) 



p"(0) = p°°(0) = |r;)H, pi"(0) = p"i(0) = 0. (34) 
An example of the master equation is presented in Section 

El 

IV. Single Photon Signal Model 

In Section Un] we saw that the master equation for the sys- 
tem G driven by a single photon field is non-Markovian, and 
the equations derived suggest the possibility of embedding 
the system and field in a larger system G. Indeed, Markovian 
embeddings were used in [9] to derive quantum trajectory 
equations for a class of non-Markovian master equations. 
In engineering and statistics, it is common practice to use 
'generating filters' driven by white noise to represent colored 
noise. 

In this section we construct a generating filter M ~ 
{I,Lo,Ho) (an open quantum system on a Hilbert space 
S)o) driven by vacuum to represent the single photon field. 
Fig. 121 Here, the ancilla parameters M = (/, Lq, Hq) are to 
be determined. This results in an extended system G defined 
on the Hilbert space S)q^ Sj (cascade, or series connection, 
[15]) driven by vacuum, with parameters given by 



G 



G<iM =iS,L + SLo, H + Ho + lm[L*SLo]), 



(35) 



from which the master equation and quantum filter equations 
(Section can be obtained. 
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Fig. 2. An ancilla system M (a two-level system) is used as a signal 
model or 'generating filter'. The ancilla system M is driven by vacuum 
(quantum white noise), and produces a single-photon output. The cascade 
circuit illustrated in this figure is equivalent to the circuit of Fig. [T] 



Specifically, we seek an ancilla system M initialized in 
a state I7) and driven by vacuum that provides a unitary 
dilation of the single photon driven system. This means that 



given our system G — {S,L,H), can we find an ancilla 
system M with state I7) and an ancilla operator R such that 
if U{t) is the unitary for G = G <i 7\/ then 



^n^Xit)] =E^„4U*it){R(g>X)U{t)] 



(36) 



Here, X{t) = U*{t)XU{t), where U{t) is the unitary for G. 
Equation ( |36] | means that the effect of the single photon field 
on the system is equivalent to the effect of the cascade of the 
ancilla on the system. The left hand side of ( |36] | is a quantum 
expectation with respect to the state \rj^) = \ri)^\£_) with the 
field in the single photon state |^), while the right hand side 
involves quantum expectation in the extended system with 
respect to the initial state |7r/0) = \'y) ^ \r]) (^ \(j>) , with the 
field in the vacuum state 10). 

In order to fulfil the requirement ( [36] l, we consider the 
time derivative of both sides of ( |36l ) and determine the 
unknown signal model parameters by comparison. The time 
derivative of the left hand side of ( |36] | is given by ( 1211 1 or 
equations (l23Tl-(l26b above, while the right hand derivative can 
be computed from the Lindblad generator Qg<im{A(E)X) for 
the extended system G — G < M, 



-E^,4U*it)iA®X)Uit)] 
E^r,4>[U*{t)gG<>M{A® X)U{t)] 



(37) 



where 



gG<M{A(E)X) = A(E)gG{X)+CL„{A)(E)X 

+L*A (g) S* [X, L] + ALo ® [L* , X] S 
+L*Lo<»{S*XS~X) (38) 

for any ancilla operator A and system operator X. 

Considering the definitions (l22l i of /Xj (X) above, we must 
find ancilla operators Qjk and (non-vanishing) weighting 
functions Wjk{t) such that 



Jk 



^^^{x) 



p.t{Qjk®x) 

Wjk{t) 



where 



MQjk ®x) = E^^^[u*{t){Qjk ® x)u{t)]. 



(39) 



(40) 



Now 



d ptiQjk^X) _ p.t{gG<iM{Q3k®X)) 



dt Wjk{t) 



Wjk{t) 

PtiQjk ®X) Wjk{t) 



(41) 



Wjk{t) Wjk{t) 

Comparing this expression with equations (l23])-(l26b we find 
that ( |39l ) and ( |36] | are satisfied if we choose the ancilla to 
be a two-level system with state 7 = |ei) (excited state), 
R = Qii = I, Qw = g-, Q oi = cr+, Qoo = n, Wii{t) = 1, 
wioit) = woi{t) = ^/w{t), woo{t) = w{t) = J^°° \^{s)\'^ds. 

The signal model ancilla system is therefore 



M^(/.«^,0), 



(42) 



and so the extended system is 



{S,L 



JMt) 



G = G<iM 



(43) 



w{t) 



Im(L*S'cr_)). 



Remark. The output state of the generating filter M can 
be understood as follows. If V{t) denotes the unitary for 
M driven by vacuum, and if the initial state is |'0(O)) = 
|ei) (g) 10), then the state lil^it)) = V(t)\ip{0)) satisfies the 
Schrodinger equation 



dm)) 



m ,„,,,, urn? . 

-a^dB [t) — - — -——aj^(T-dt 



VW) 



2 w{t) 



mt)) 

(44) 

(since di?(t)|(/)) = 0). It is an elementary calculation to see 
that this has the exact solution 



m)) - y^i)|ei ® 0) + |eo ® sr*(e)0); 



(45) 



cf. ([8]l. Since ^ is square integrable, |^(t)) approaches the 
state \eo (E) £,) as t —>■ oo. Therefore the state of the field 
produced by the generating filter M approaches the single 
photon state |^) asymptotically. 

V. Single Photon Filter 

The quantum filter Trt{X) = E^^[X(i)|y(s), < s < t] 
for the system G driven by a single photon field may 
now be obtained from the quantum filter TCtiA (g) X) = 
E^r,4U*{t){A(g) X)U{t)\I (g)Y{s),0 < s < t] for the 
extended system G = G < M driven by vacuum [8], with 
the extended system parameters given by ( |43] l (see the 
Appendix). Using the definition of conditional expectation 
(see [8] and the Appendix), it follows that 



I(g)Trt{X)^nt{I(g)X). 



If we define 



jk. 



I®<{X) 



ntiQjk<^X) 

Wjk{t) 



(46) 



(47) 



and the innovations process W{t) is a Wiener process with 
respect to the single photon state and is defined by 



we obtain the system of equations 

dnl^X) 
= {nl\giX))+7T^\S*[X,L])Cit) 

+7Tl°i[L*,X]S)m + ^^'iS*XS~X)\m\')dt 
Ht^IHXL + L*X) + 7r°\S*X)C{t) + 7:l°{XS)m 
-7Tl\X)ICt{X))dW{t), (48) 

dnl'iX) = {nl"{g{Xj)+7r'i°{S*[X,L])C{t))dt 
+{7rl°{XL + L*X) + Tr°^{S*X)Cit) 
-TTl°{X)ICt{X))dW{t), (49) 



dn^iX) 



T:'>"{g{X))dt + «'(XL + L*X) 
^7T°''{X)ICt{X))dWit). (50) 



Here, 



ICtiX) = 7:l\L + L*) + 7:riS)m + 7rr(^*)r W (51) 



dW{t) ^dY{t)~]Ct{X)dt. 



(52) 



We have 'n^^{X) = nl^iX*)*, and the initial conditions are 
^l\X) = nf>{X) = (r?,X,y), nl^{X) = ^§i(^) = 0. 

Consequently the conditional expectation for the system 
driven by the single photon field is given by 

nt{X) = K^^[X{t)\Y{s), < s < t] = n]\X), (53) 

and so the required quantum filter is given by the system of 
coupled equations (|48]l-(l50ll. 

Equations for the conditional density operators may easily 
be derived. Finally, we remark that although the master 
equations for systems driven by a single photon field have 
been obtained by other means in [13], to the best of our 
knowledge the quantum filtering (trajectory) equations for 
such systems have not been derived before. 

VI. Example 

When the system is a two-level system or qubit, the filter- 
ing equations reduce to a finite set of stochastic differential 
equations. In this case we have S^s = C^. The system 
is specified by the parameters S ^ I, L = ^/ku-, and 
11 = uja^. Here k > is a scalar parameter. 

We begin with the master equations (I30l) - (l33b . and write 

^" (54) 



„oi 



^11 



-(J + x^V. + y'^V. + z'^V,) 

(a;"V, + y"V, + z"V,)-p 



(j + X^V^+y^Vy + Z^V 



(55) 
(56) 



Note that x^^ , -ip^ , z"" and x^^, y^^, z^^ are real, while x^^ , 
2/°^, z°^ may be complex. Also note, for example, p^^ia^) = 
X™, p°^((Tj;) = x^^* , etc. Then we obtain nine coupled 
equations for the nine coefficients: 



X 



00 



-2c.y™-fx°°, 



■00 



„00 



y"" = lujx^" - -y 

iOO 



00 



Z^" = -K.(l + Z^^), 



X 



01 



x"i-2wyOi-y;^^(t)*z 



*„00 



y"i = 2..x"i-^z/i-*V^C(t) 



*^oo 



iOl 



■,11 



KzOi - V^xO"e(i)* + i%/^y°°e(i)*, 



'2" 

11 



2a;x 



2wy" + V^z°iC(i) + 4^z''^*i{tY 



i" = -K - Kz" - V^x°i^(i) - iy^y^'^^it) 

For the quantum filter (l48])-(l50b. we use a slightly more 
general representation for fp^ given by: 

1 



P' = T^ic' I + S:' (Tx + y^'^ay + 



?J'^. 



for j, fc = 0, 1. Since {p^^,I) = 1 (i.e., p^^ is a normalized 
conditional density operator), we always have that c}^ ^ 1 
at all times. However, unlike the master equation, this will 
not be so for cP^, c^", c"", as these coefficients will evolve in 
time. This is the reason we need to consider the more general 
representation for p^'^. The quantum filter for the two-level 
system is given by the finite set of coupled equations 



rfgOo^r/T^oo 



(T^jii + l£Oi^(^) ^ h'^*^{ty)c°°)dW{t), 



di"° - (-2c^y"° - ^x°")dt + ( V^c 



,00 



(V^x" + ic°ieW + ^c°i*CW*)i™)rfW^(i), 






dx" = (-|.t" - 2c^y" + V^z"^C(i) + V^z°^*atr)dt 



+ (v^+x"^*atr+i°^^{t) 



(^T^ii + ic°icw + ic"i*cW*)-i")dw^W, 



#' 



2a;a;' 



-y" + zV^2°ie(i) - iV^z°'*eW*)di 



f'*m* + rm 



(y^jii + ^c"'m + lc"'*m*)f'ywit), 



dz^^ = {-K - Kz" - ^^x°l^(t) - ^^/^y"'^(i) 



The innovations process is given by 

dW{t) ^ dY{t)- {^x^\t)+c''\t)m+c^\t)C{t))dt. 

(57) 



dc"i= f^^x"l+c"°CW* 



(y^jii + i£Oi^(i) + ic"i*^(i)*)c°i)dTy(f), 






dy" 



(y^^ii + lc°i^(t) + ic°i*eW*)i°')dVr(t), 



(2^^oi_|yOi_^^^(^)*^oo)^^ 



dz" 



rC(t)' 



(T^^ii + icOi^(t) + \c'^'*m*)f')dW{t), 



■ (-Kz"i - y^x™C(t)* + iy/^y°"£_{ty)dt 



(7^:^11 + icOi^(t) + ^c''^*atr)z°')dW{t), 



VII. Discussion and Conclusion 

In this paper we have derived the master equation and 
quantum filter for a class of open quantum systems that are 
coupled to single photon fields. The paper has focused on 
the case of quadrature measurements Y{t) given by (fT4l i. 
(flSl l. However, the methodology also applies to the case of 
counting measurements, corresponding to a photodetector in 
place of the homodyne detector in Figure [To 

The single photon filter consist of coupled equations 
that determine the evolution of the conditional state of the 
system under continuous (weak) measurement performed on 
the output field, in contrast to the familiar single filtering 
equation for open Markov quantum systems that are coupled 
to coherent boson fields. This coupled equations structure 
of the master and filter equations is a reflection of the non- 
Markov nature of systems coupled to single photon fields. 
Indeed, a key feature of our approach is the embedding 
of the system into a larger extended system, a technique 
often employed in the analysis of non-Markov systems, 
providing an elegant framework within which to study the the 
dynamics, both unconditional and conditional, of the system. 
We expect that the basic approach taken in this paper can be 
adapted to study quantum systems that are coupled to other 
types of highly non-classical boson fields. 

'^Of course, homodyne detection is based on a photon counting system, 
e.g. [12]. 



Appendix 

The quantum filter for a system G = {S, L, H) (driven 
by vacuum) for the quadrature output field observable Y{t) 
(given by (fT4l i) is 

d7:t{X) ^ 7rt{gG{X))dt (58) 

+{TrtiXL + L*X) - TTt{L + L*)ntiX))dW{t), 

where W{t), a Wiener process called the innovations pro- 
cess, is given by dW(t) = dY{t) - TTt{L + L*)dt). By the 
spectral theorem it follows that the quantum filter is equiva- 
lent to a classical system driven by the measurement record, 
[7], [8]. Perhaps the simplest way to derive the quantum filter 
is to use the conditional characteristic function, [5], [7], [25], 
which we now briefly summarize. 

We begin with a short discussion of quantum conditional 
expectation. The measurement signal Y{s), Q < s < t, 
is a collection of commuting self-adjoint operators. These 
operators form a subspace ^ in the space of operators, 
and the quantum conditional expectation ¥,,-i^[X{t)\Y{s), 
< s < i] is the orthogonal projection of the system 
operator X{t) at time t (since the field serves as a probe, 
we have the commutation relation [X{t),Y{s)] ~ for all 
< s < i (non-demolition), and so the conditional expec- 
tation is well-defined). The orthogonal projection property 
corresponds to least squares estimation, and leads to the 
following characterization ( [TtI i mentioned in Section [III We 
will use this characterization in the following form. Define, 
for any function g, 

cg{t) = cxp{ / g{s)dY{s) -If \g{s)\'ds}, 
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1 

2_ 
and note dcg{t) = g(t)cg{t)dY{t). Then we require 

E^4X {t)Cg{t)] = Er,4,[X {t)Cg{t)l 

for all functions g. 

Now suppose that dX{t) has the form 



(59) 



(60) 



(61) 



dX{t) =a{t)dt + l3{t)dY{t), 

where a and /3 are to be determined from the relation ( l60b . 
Now using the QSDE ^ for X{t) = jt{X), we have 

iE„4X{t)Cgit)] 

^En4cg{t)7rti-i[X,H]+gL{X)) 

+g{t)cg{t)TTt{XL + L*X)] (62) 

(here we have used property (l60ll). Similarly, using dMT l we 
have 

iE,4Xit)Cgit)] 

= E,j4cg{t){a{t) + nt{L + L*)p{t)) 
+g{t)cg{t){m + TTt{L + L*)X{t)] (63) 

Now equating the RHS of (|62] | and ( |63] | and using the fact 
that g is arbitrary we find that 

ait) = 7rti~i[X,H]+gLiX))-TTtiL + L*)l3{t) 

/3(t) - 7rt{XL + L*X)-TTtiL + L*)TTt{X) 

The quantum filter (ISST i now follows from this. 
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